Selected topics in homogenization of transport processes in 
historical masonry structures 

Jan Sykora a , Jan Zeman a ' b , Michal Sejnoha a '* 

a Department of Mechanics, Faculty of Civil Engineering, Czech Technical University in Prague, 
Thdkurova 7, 166 29 Prague 6, Czech Republic 
h Centre of Excellence I T4 Innovations, VSB-TU Ostrava, 17. listopadu 15/2172 708 33 Ostrava-Poruba, 

Czech Republic 



Abstract 

The paper reviews several topics associated with the homogenization of transport pro- 
cessed in historical masonry structures. Since these often experience an irregular or ran- 
dom pattern, we open the subject by summarizing essential steps in the formulation of a 
suitable computational model in the form of Statistically Equivalent Periodic Unit Cell 
(SEPUC). Accepting SEPUC as a reliable representative volume element is supported by 
application of the Fast Fourier Transform to both the SEPUC and large binary sample of 
real masonry in search for effective thermal conductivities limited here to a steady state 
heat conduction problem. Fully coupled non-stationary heat and moisture transport is 
addressed next in the framework of two-scale first-order homogenization approach with 
emphases on the application of boundary and initial conditions on the meso-scale. 

Keywords: Random masonry, Binary images, Statistically equivalent periodic unit cell, 
Computational homogenization, Fast Fourier Transform, Coupled heat and moisture 
transport 



1. Introduction 

Advanced computational simulations of historic structures are becoming increasingly 
common in the assessment of their existing state and in planning of reconstruction mea- 
sures [10]. In this context, particular attention needs to be paid to variations of tempera- 
ture and moisture fields, whose contribution to structural damage usually far exceeds the 
effects of mechanical loadings, see e.g. [19, 31, 2] for concrete case studies. Taking into 
account highly heterogeneous character of historical constructions, simulations of these 
phenomena often necessitates the deployment of multi-scale strategies, developed for me- 
chanical [1, 11, 12, 28] and transport processes [24, 25] in masonry structures. Successful 
engineering application of this modeling approach to the assessment of Charles bridge in 
Prague is described in our work [33]. There, three-dimensional multi-physics analysis of 
the bridge body was executed, with material parameters at the macroscopic structural 
scale determined from meso-scale simulations. Among other things, results of the study 
highlighted the need for fully coupled macro-meso simulations of transport processes. 
In this paper, we complement these results by a more detailed analysis of two aspects of 
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multi-scale simulations, namely the meso-scale representation for irregular masonry struc- 
tures and the introduction of proper boundary and initial conditions in the macro-meso 
scale transition. 

The masonry texture at meso-scale is to be incorporated in the form of the Repre- 
sentative Volume Element (RVE), a statistically representative sample of the analyzed 
material [4]. Specification of the RVE is particularly simple in case of a regular ma- 
sonry, for which it reduces to a Periodic Unit Cell (PUC) associated to a given type of 
the bond [1]. However, historical masonry structures are typical of irregular or random 
textures, which renders the determination of PUC unambiguous. A convenient approach 
to overcome this difficulty was introduced by Povirk [18], who suggested to replace the 
original complex meso-structure with an idealized PUC, with parameters determined by 
matching spatial statistics of original and simplified representation. By combining his 
ideas with related works on microstructure reconstruction, e.g. [32], the concept of Sta- 
tistically Equivalent Periodic Unit cell (SEPUC) was later successfully applied to, e.g. 
analysis of fiber-reinforced [35, 30] and woven composites [36, 27], high-density poly- 
disperse particle packings [9] or micro-heterogeneous steels [21]. In Section 2, we employ 
this procedure to an irregular masonry wall of Charles Bridge in Prague, characterized 
by a digital photograph. 

Generating a suitable RVE of a heterogeneous material is just the first step towards a 
reliable prediction of material as well as structural response of masonry. Focusing on the 
description of transport processes in heterogeneous media the reader is advised to study 
the work by Ozdemir et al. [16] and Larsson et al. [8]. Through the application of consis- 
tent variational formulation the authors in [8] suggested the macroscopic response to be 
dependent on the actual size of mesoscopic RVE, providing the transient conditions are 
assumed on both the macro- and meso-scale. Being aware of the need for a fully coupled 
multi-scale analysis of simultaneous heat and moisture transport in masonry structures 
performed recently in [23] , it becomes clear that additional sub-stepping of a given macro- 
scopic time step on the meso-scale may considerably slow down the computational process. 
To show that in some cases this step might be avoided by running the meso-scale analysis 
under steady state conditions [16] even for a finite size RVE may thus prove useful in 
keeping the computational cost relatively low. This issue is addressed in the second part 
of this paper, Section 3, with particular attention dedicated to the influence of loading 
and initial conditions imposed on meso-scale. 

2. Image-based geometrical modeling 

This section deals with the first aspect of multi-scale simulations announced above, 
namely with the realistic representation of irregular masonry structure relying on image- 
based data. To this purpose, in Section 2.1 we introduce the model of an idealized 
mesostructure and the procedure to determine its parameters. Section 2.2 briefly reviews 
the numerical scheme employed to determine local fields operating directly on mesostruc- 
tural images. These two tools are combined together in Section 2.3 to assess the SEPUC 
quality in view of the distribution of heat fluxes in disordered masonry under steady-state 
conditions. 

2.1. Strategy of SEPUC determination 

The key step in the SEPUC definition is a proper choice of the spatial statistics to 
characterize the dominant features of the heterogeneous material under study. With the 
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focus on masonry structures, we limit our attention to the mortar two-point probabil- 
ity function, see e.g. [26] for a general overview and examples of alternative statistical 
descriptors. 

To introduce the subject, consider a masonry sample Q composed of mortar and blocks, 
and denote the characteristic function of the domain occupied by mortar as x( x )- Then, 
the two-point probability function £2 states the probability that two points x and 
randomly thrown into a medium, will both be found in the mortar phase: 

S 2 (x,y)=F(x(x) X (y) = l). (1) 

For the case of statistically homogeneous and ergodic media, two-point probability func- 
tion depends on (x — y) only and 62(0) = 0, where (j) is the mortar volume fraction. 
Moreover, it can be obtained from the relation 

= y^F- 1 (7 (x)7W) , (2) 

that can be efficiently evaluated using the Fast Fourier Transform (FFT) techniques for 
image-based microstructures [3, 26]. In Eq. (2), T (•) and J 7-1 (•) designate the forward 
and inverse Fourier Transform operators, respectively, • denotes the complex conjugate 
and is the area of Vt. 

Once the original structure has been quantified by a suitable statistical descriptor, a 
proper parametrization of the idealized cell geometry needs to be introduced, expressed 
here by a parameter vector p. Its optimal value then follows from the minimization of 
least square error 

E W ~ TcT7 f B x ) " s ^P)f dx > ( 3 ) 

expressed as the difference between the target statistical descriptor £2(03) related to the 
original microstructure and the unit cell associated with p, integrated over the unit cell 
domain Q n C 

A closer inspection reveals that the objective function (3) is non-convex, multimodal 
and discontinuous due to the effect of limited bitmap resolution. Based on our previous 
experience [5, 13], a global stochastic optimization algorithm, relying on the combination 
of real- valued genetic algorithms and the Simulated Annealing method, is employed to 
solve this optimization problem. 

2.2. Homogenization scheme 

The distribution of local fields within Q n follows from the solution of periodic unit cell 
problem [14] 

VxV# = 0, V-q = 0, q = -XV9 in ft D , (4) 

in which q stands for the thermal flux vector, A is the second-order tensor of material 
conductivity and V# denotes the ^-periodic field of temperature gradient satisfying 

/ V6(x)dx = V6, (5) 

where VB is the macroscopic temperature gradient prescribed over Q n , see Section 3 
below for additional details. It is well-known that the solution to the unit cell satisfies 
the Lippman-Schwinger equation [14] 



V6(x)+ [ T(x-y)SX(y)Ve(y)dy = Ve, 



(6) 
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where 5X = X — X^I, A is the conductivity of an auxiliary isotropic reference medium 
and the second-order operator T is related to the Green function of the problem (4) with 
X(x) = A(°) J. It admits a compact closed-form expression in the Fourier space [14] 

fo forfc = 0, 

T(T)(k) = {k®k ^ . (7) 

— — — otherwise, 

K k - k 

so that its action can be efficiently evaluated by the FFT algorithm. This observation 
resulted in an iterative scheme due to Moulinec and Suquet [15], applicable to arbitrary 
digitized media. 

In our case, we adopt an accelerated version of the original algorithm due to Zeman 
et al. [34]. Since the sample is discretized by a regular N\ x N 2 bitmap, it is convenient 
to project the integral equation (6) onto the space of trigonometric polynomials, e.g. [20]. 
This yields the linear system in the form 

(I + B)V0 d = V0 d , (8) 

where the 2NiN 2 vector V#d stores the unknown discrete values of temperature gradient 
at pixels, VOd is the corresponding vector of the overall temperature gradient and matrix 
B is expressed as 

(F- 1 \ (T n T 12 \ (F 0\ fSX n SX 12 
^ V F- 1 ) \T 21 T 22 ) \0 FJ \SX 21 SX 22 

Here, the matrices F and F~ x implement the forward and the inverse discrete Fourier 
transform and, e.g. 8X\ 2 is a diagonal (NiN 2 ) x (NiN 2 ) matrix storing the corresponding 
component of the conductivity matrix at individual pixels, see [34] for more details. The 
system (8) is solved using standard conjugate gradient algorithm. 

Upon convergence, the distribution of the local heat flux q due to VB is determined 
from the solution V#d by Eq. (4) 3 and the global heat flux is computed as 

Q = w~\ f q ^ dx - ( 10 ) 

I I— 1 1 J £2\J 

This allows us to determine the 2x2 homogenized conductivity matrix A hom via the 
solution of two successive steady state heat conduction problems. To that end, the periodic 
unit cell is loaded, in turn, by each of the two components of VB equal to unity, while 
the remaining one vanishes. The corresponding volume flux averages Q then provide 
individual columns of A hom . 

2.3. Example 

Principles of the introduced methodology are illustrated by the analysis of sandstone 
facing masonry wall of Charles Bridge in Prague appearing Figure 2.3. As the first step, 
the original color image, Figure 2.3(a), was thresholded to the binary representation, 
manually adjusted to remove image processing artifacts and a rectangular domain Q was 
selected for further analysis, Figure 2.3(b). 

Our aim is to replace this complex meso-structure with its idealized representation 
in terms of a two-block-layer SEPUC, described by twelve parameters assembled in the 
vector 

P= (6,61,62, A,/li,/l 2 ,ti,t 2 ,t3,t4,t 5 ,t 6 ) , (11) 
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(a) 



(b) 



Figure 1: (a) Example of a masonry wall of Charles Bridge in Prague and (b) meso-structural window Q 
provided in the form of 1,600 x 735 black- and- wide bitmap 

see Figure 2(a). These parameters are adjusted to minimize the discrepancy between 
the target two-point probability function, Figure 3(a), and the one corresponding to 
the SEPUC, Figure 3(b). As visible, the resulting SEPUC-based representation Q n , 
Figure 2(b), captures the dominant features of the original microstructures, such as the 
overall volume fraction equal to 7.5%, or the average thickness of joints, corresponding 
to the spread of the peak at x = in the horizontal and vertical directions. However, 
a significant amount of fine-scale features of the two-point probability functions, visible 
in the original data, are filtered out or poorly reproduced by the adopted, more regular, 
representation. 

Of course, the impact of such approximation needs to be quantified a-posteriori, from 
the point of view of the particular application. In our case, such a step is performed on 
the basis of homogenized conductivities and distribution of heat fluxes as determined by 
the FFT scheme. The distribution of the magnitudes of local fields due to macroscopic 
temperature gradient dO/dx 2 appears in Figure 4, for phase conductivities set according 
to Table 2. We observe that even though the local fields within the SEPUC are certainly 
more regular than in the original mesostructure, the extreme values as well as the average 
distributions are reproduced surprisingly well. 

This claim is further supported by very close match between the homogenized con- 
ductivities of SEPUC and the original media, Table 1, demonstrating that the SEPUC is 
capable of reproducing almost perfect isotropy of the original sample. Finally, for the sake 
of comparison, we also present homogenized conductivities of unit cells used in the fol- 
lowing section, cf. Figure 5. The influence of the degree of heterogeneity on the resulting 
predictions is evident. 




horn \ horn 



21 ^22 



Original mesostructure 0.2612 



0.0000 0.2622 
0.0000 0.2618 



SEPUC 0.2616 



Regular masonry 0.2856 
Irregular masonry 0.2932 



0.0000 0.2967 
-0.0015 0.2960 



Table 1: Homogenized thermal conductivities in [Wm X K x ] 
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Figure 2: Determination of SEPUC for irregular masonry; (a) adopted parametrization and (b) optimal 
mesostructure with dimensions in [mm] 
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Figure 3: Two-point probability function describing (a) target microstructure and (b) SEPUC (1 pixel 
corresponds to ~ 2.6 mm) 
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Figure 4: Distribution of heat flux magnitudes \\q\\ (in Wm 1 ) due to due to macroscopic temperature 
gradient dO/dx2 = 1 Km -1 in (a) original domain and (b) SEPUC 
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3. Multi-scale homogenization of coupled heat and moisture transport 



In this section, we continue with the selected example of Charles Bridge and extend 
the previous study to the multi-scale modeling of coupled nonlinear transient heat and 
moisture transport in masonry structures. While still adopting the first order homog- 
enization approach (linear variation of macroscopic temperature and moisture fields is 
assumed) we choose, unlike the previous section, the finite element method (FEM) to 
solve the resulting system of partial differential equations. In the present study, these 
arise from the application of a nonlinear diffusion model proposed by Kiinzel in [6, 7]. 
The model is described by the energy balance equation 



dH d6 
~d6di 



= V T [AV#] + h v V T [5 p V{(pp sat (6)}} 



(12) 



and by the mass conservation equation 



dw dip 
dip dt 



V T [D 9 Vip] + V T [<5 p V{ mat (0)}] 



(13) 



where temperature 9 and relative humidity cp are the two state variables. The model 
requires specifying the enthalpy of the moist building material iif , the water content of 
the building material the coefficient of thermal conductivity A, the liquid conduction 
coefficient D ip) the water vapor permeability 8 P) the evaporation enthalpy of water h v and 
the water vapor saturation pressure p sat . Details regarding the functional dependence of 
the model parameters on temperature and relative humidity are available, e.g. in [6, 22]. 



b' 




in = 14[°C] w 

< \ 'X I ^ 1 -X I 



-k = 1.92 [mi- 



Mac ro-sc ale 



1, = 0.32 [m] 

Meso-scale 



Figure 5: Scheme of coupled multi-scale framework 

In [22] the authors noticed a significant dependence of the homogenized macroscopic 
properties on the applied temperature and relative humidity gradients. This finding 
promoted the solution of full scale analysis of masonry structures in a fully coupled multi- 
scale homogenization framework. To that end, a nested FE 2 scheme, graphically presented 
in Figure 5 (see also [16, 8, 23] for more details), appears as a suitable method of attack. 
In such a case we expect the homogenized macro-scale fields to be found from the solution 
of a certain sub-scale (meso-scale) problem performed on an RVE loaded by the prescribed 
constant temperatures and moisture gradients. 
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As already mentioned in the introductory part, Larson et al. [8] have shown that 
assuming a non-stationary response also on meso-scale generates a non-local term in the 
homogenized macroscopic equations, which renders the macroscopic response dependent 
on the RVE size. There is no dispute that performing the non-stationary analysis on both 
scales may considerably increase the computational cost making the numerical analysis 
prohibitively expensive even if exploiting parallelization. To reconcile these issues for 
a typical masonry material and RVE sizes is thus crucial for the success of analysis of 
full scale tree-dimensional models of historical masonry structure such as Charles Bridge. 
This will be the main topic of subsequent paragraphs addressing also the effect of loading, 
boundary and initial conditions. 

Since irrelevant from the computational point of view we consider, henceforth, only 
the regular (PUC) and irregular (SEPUC) bonding of masonry displayed in Figure 5. 

3. 1 . Theoretical formulation 

To begin, we adopt a variationally consistent homogenization outlined in detail in [8] 
starting with the assumption that a local field a can be replaced by a spatially homogenized 
one (a) such that 

/ adn « / (a) u dQ= [ ( -J— [ adft D )dft, (14) 

/ adr ~ L^ ndT = LiwiL adra ) dr ' (15) 

where Q n and T n represent the internal and boundary parts of PUC. In what follows, 
for the sake of lucidity, we shall treat only the energy balance equation (12) which upon 
employing Eqs. (14) and (15) becomes 

L «^) Aa + L ( { ™ }T [ { AV " } + K ™}] ) a in+ 

+ / ({vse} T [h v {s pPs , t v^}]) n dn- [ <^) D dr = o. (ie) 

JQ 1 JY q e 

In the spirit of the first order homogenization, it is assumed that the macroscopic 
temperature and relative humidity vary only linearly over PUC. This can be achieved 
by loading its boundary by the prescribed temperature O hom and relative humidity $ hom 
derived from the uniform macroscopic temperature VB and relative humidity V$ gradi- 
ents. In such a case, the local temperature and relative humidity inside PUC admit the 
following decomposition 

0(x) = G(X°) + {VO} T {x-X°} + 9*(x) = G hom (x) + 9*(x), (17) 
<p(x) = <5>(X ) + {V<5>} T {x-X°} + ^(x) =$ hom (^) + (^*(^), (18) 

where 0*(x) and (p*(x) are the fluctuations of local fields superimposed onto linearly 
varying quantities Q hom (x) and $ hom (cc) . The temperature B(X°) and the moisture 
at the reference point X° are introduced to link the local fields to their macroscopic 
counterparts. For convenience, the PUC is typically centered at X°. Henceforth, the 
local fluctuations will be demanded to be periodic, i.e. the same values are enforced 
on the opposite sides of a rectangular PUC. This ensures the scale transition condition, 
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see e.g. [16], which arises from upon averaging the micro-temperature gradient over the 
volume |0| of PUC 

({vo(x)})n = -r±- [ {ve(x)}dsi(x) = {ve} + ^- f {ve*(x)}<m(x). (19) 

J Sin J Sin 

Hence we demand the contribution of fluctuation fields to disappear upon volume aver- 
aging 

<{W(z)})o = ^ / {Vd*(x)}dQ(x) = -±- [ P(x){u(x)}dT{x) = 0, (20) 

where 1/ stores the components of the outward unit normal to the boundary of PUC being 
directed in opposite directions on opposite sides of the PUC. 

Next, substituting Eq. (17) into Eq. (16) and collecting the terms corresponding to 
56 hom and 86* splits the original problem (12) into the homogenized (macro-scale) prob- 
lem and local sub-scale (meso-scale) problem. Since details on the derivation of equa- 
tions driving the solution on individual scales have already been given in our proceeding 
works [25, 23], we present only the result pertinent to the macro-scale 

-/> e >lfs( e+veT {*-* o })) D d!i - 

V v ' 

- £ ( { sveV {x - x»}f A (e + ve> - *■})) d n - 



J 



(c ee +c ee ^ 



dt 



{sve} T 

n 



dft- 

□ 



- / ({5ve} T [h v {5 pPsat v^}]) n dn+ [ ({5e} T ^> n dr = o. (21) 

K e(p r v g t 

to identify the solution dependence on the actual size of PUC through the second term 
in the integral (21). 

An analogous approach can be applied also to the moisture transport equation (13) 
to arrive, after classical finite element discretization, into a discretized system of coupled 
macroscopic heat and moisture equations 

K^ + K^ + (C„ + 2C;, + C;;)^ = </ ext , (22) 

K^ + K w r^ + (C W + 2C W + C^)^ = 9ext , (23) 

which have to be properly integrated in the time domain adopting for example the 
Crank-Nicolson integration scheme. Details on the numerical implementation are avail- 
able in [22]. The homogenized matrices in Eqs. (22) and (23) follow directly from the 
meso-scale solution for a given macroscopic time increment. Because of a strong non- 
linearity arising on both scales the two problems (macro-meso) must be solved iteratively 
by the Newton-Raphson method, see [16, 8, 23] for further reference. 
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3.2. Numerical examples 



A fully coupled multi-scale analysis of a two-dimensional segment of Charles Bridge 
in Prague has been carried out in [23] with emphases on parallel computing. This study 
will be extended herein by quantifying the influence of the boundary, loading and initial 
conditions when moving down from the macro-scale to the meso-scale. 



3.2.1. Boundary conditions 
To begin, recall Eq. (20) 

<{W (*)})□ = J-r / P(x){v(x)} dT(x) = 0, (24) 

being satisfied providing either the fluctuation part of the temperature field equals zero 
or the periodic boundary conditions, i.e. the same values of #* on opposite sides of a 
rectangular PUC, are enforced on r n . 



• -^z^ 

"a- ya 
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X 

> D (master) 
O ( 


) o 



Figure 6: Periodic boundary conditions 



Such conditions are easy to impose if searching the solution in terms of the fluctuation 
part of the temperature or moisture fields. The macroscopic constant gradients {VB} 
and {V$}, see Figure 7(a), are then directly used to load the unit cell [25]. However, 
this is mostly not possible with the application of commercial codes. In such a case, the 
constant gradients are introduced by enforcing a linear variation of the homogeneous part 
of local fields 6 hom (cc), $ hom (£c). This is achieved by prescribing directly the Dirichlet 
boundary conditions along all edges of PUC, i.e. 9(x) = Q hoin (x),(p(x) = $ hom (cc) (see 
Figure 7(b)) are specified on the boundary r n of PUC. The periodic boundary conditions 
are then prescribed with the help of multi-point constraints. 

In doing so, observe in Figure 6 that for a pair of points (e.g. A - master and a - slave) 
located on the opposite sides of PUC the following relations hold: 

eA=(j?)vA + P A + 9(x ), Va=(^) L +(^)y« + e: + ®(x ). (25) 

Taking into account the fact that the fluctuation field 9* satisfies the periodicity condition 

K = Pa, (26) 
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and subtracting corresponding terms on the opposite edges, we finally obtain (compare 
with [16]) 

L = 8a — 9A = 92 — 9i, 
(^)# = 6 b -6 B = e z -6 u (27) 

where 1? 9 2 and # 3 are the temperatures at the control points 1,2 and 3 seen in Figure 6. 
The same conditions apply to the relative humidity tp as well. 

3.2.2. Loading and initial conditions 

An illustrative example presented here considers a finite element mesh consisting of 108 
macro-elements each representing a single meso-problem with assigned periodic boundary 
conditions. Its geometry together with the specific macroscopic loading conditions are 
shown in Figure 5. 

The initial temperature 9- m = 14 [°C] and the moisture ip- m = 0.5 [-] were assigned to 
the whole domain. The following boundary and loading conditions were imposed: the left 
boundary of the domain was submitted to exterior loading conditions ext = 5 [°C] and 
9?ext = 0.5 [-], while the opposite side was submitted to interior loading conditions 9- int = 24 
[°C] and (pint = 0.8 [-]. Zero flux boundary conditions were assumed for horizontal edges. 



parameter 






brick 


mortar 


Wf 


[kgm" 


- 3 ] 


free water saturation 


229.30 


160.00 




[kgm" 


- 3 ] 


water content at ip = 0.8 [-] 


141.68 


22.72 


A 


[WHI- 


-i K -i] 


thermal conductivity 


0.25 


0.45 




CH 




thermal conductivity supplement 


10 
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Ps 


[kgm" 


- 3 ] 


bulk density 


1690 


1670 


fl 


["] 




water vapor diffusion resistance 


16.80 


9.63 


A 


[kgm" 


-2 s -0.5] 


water absorption coefficient 


0.51 


0.82 




[jk g - 




specific heat capacity 


840 


1000 



Table 2: Material parameters of individual phases 



The material parameters in Table 2 were obtained from a set of experimental measure- 
ments providing the hygric and thermal properties of mortars and bricks/stones, which 
have been used in the reconstruction works of historical buildings in the Czech Republic, 
see [25, 17]. 

With reference to Section 3.2.1 we consider two types of loading conditions. Fig- 
ure 7(a) assumes the loading in terms of macroscopic temperature and moisture gradients 
to provide local fluctuations upon submitting decompositions (17) and (18) into Hill's 
averaging condition (consider for simplicity a steady-state heat conduction problem only) 

({SVd} T {q}) = 0, (28) 

to get 

({5V6*} T [A] {V0*}> = -({5V#*} T [A] {V6}). (29) 

The periodic boundary conditions are then easily introduced by prescribing the same 
code numbers, in the FEM language, to the homologous nodes on the opposite faces of 



11 



the rectangular unit cell. Figure 7(b) then presents fully prescribed Dirichlet boundary 
conditions applicable with commercial codes when searching the solution directly in terms 
of the local fields 0, (p. The periodic boundary conditions are then enforced indirectly 
using Eq. (27). 



e hom ( X ) = B(x°) + {ve} T {x-x u } 




(a) (b) (c) 



Figure 7: Different loading schemes of PUC on meso-scale: (a) prescribed macroscopic gradients (loading 
conditions I, labeled as NST 1 ), (b) prescribed macroscopic temperatures and relative humidities (loading 
conditions II, labeled as NST 2 ); (c) Assumed initial conditions 

Apart from loading conditions, the two types of initial conditions seen in Figure 7(c) 
are also examined. In particular, either a linear variation of the homogeneous part of local 
fields 6[ n (x) = 6 ° m (a0, <Pm(&) = $ ° m ( x ) is assumed or the macroscopic temperature 
^in(^) = 6(X°) and moisture (p[ n (x) = $(X°) calculated at the reference point X° at 
the end of the current macroscopic time step are assigned to the whole PUC. Note that 
the former one arise naturally from Eqs. (17) and (18) when setting the fluctuation fields 
equal to zero at time t = 0. 

Meso-scale analysis 

In the first example attention is dedicated to the solution of a coupled transient heat 
and moisture problem solely on the meso-scale. First, the regular periodic unit cell, 
PUC in Figure 5, was loaded by the highest temperature and moisture gradients ob- 
tained in the course of actual multi-scale analysis discussed later in this section. The 
following numerical values were considered; initial conditions (B(X°) = 17 [°C], = 
0.15,0.55,0.85 [-])and boundary conditions ({V6} = {V x e,V y 6} = {5.0,1.0} pCm" 1 ], 
{V$} = {V x $, V y $} = {0.1, 0.05} [m- 1 ]). 

Figure 8 shows evolution of local temperature and moisture fields along the horizontal 
centerline of PUC as a function of time derived from two types of loading conditions. Pe- 
riodic boundary conditions and linearly varying initial conditions with = 0.55 [— ] 
were considered. It is shown that for the adopted extreme gradients the steady state solu- 
tion, plotted as solid lines, is reached in about 15 hours irrespective of loading conditions. 

Further to this subject, we also suggest invariance of the solution to the assumed initial 
conditions as depicted in Figure 9. This becomes evident once realizing a nonlinearity 
of Eqs. (12) and (14) taken into account through the application of Newton-Raphson 
iteration scheme. Clearly, the initial difference in the solution error attributed to the 
difference in initial conditions is wiped out already in the first load (time) increment 
upon arriving at equilibrium. 

The influence of the degree of material nonlinearity of the present constitutive model 
is partially seen in Figure 11 showing evolution of the selected effective material properties 
as a function of time for three values of initial macroscopic moisture Functional 
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Figure 8: Evolution of local fields (PUC): (a,b) temperature, (c,d) moisture; (a,c) Loading conditions I 
(Figure 7(a)), (b,d) Loading conditions II (Figure 7(b)) 
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Figure 9: Influence of initial conditions (constant or linearly varying homogeneous part of local fields); 
Evolution of effective moisture terms as a function of time: (a) effective conductivity term, (b) effective 
storage term 



dependence of some material parameters on moisture is plotted in Figure 10 for illus- 
tration. The complete set is available, e.g. in [22, 29]. It is evident that for effective 
properties the time to reach the steady state solution is considerably shorter than for lo- 
cal temperatures and humidities, recall Figure 8. It is even more important to realize that 
the difference between effective steady state parameters and effective parameters withing 
a transient regime is essentially negligible. This can be attributed to a relatively small 
difference between local fields pertinent to steady state and various stages of transient 
solutions as evident in Figure 8. In other words, the material nonlinearity observed in 
Figure 10 does not play in this case a significant role, thus promoting the steady state 
solution from the macroscopic point of view as sufficiently accurate. 
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Figure 10: (a) Variation of water content as a function of relative humidity, (b) Variation of phase 
moisture storage function as a function of relative humidity, (c) Variation of slope of saturation pressure 
as a function of temperature 
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Figure 11: Evolution of homogenized material parameters as function of time (PUC): (a) ^ ^, (b) 
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Finally, to address the influence of the degree of material heterogeneity we repeated 
the same study employing the irregular periodic unit cell (SEPUC in Figure 5). Fig- 
ures 12(a) (b) show evolution of relative humidity along the SEPUC centerline at various 
times clearly identifying the material boundaries. While the evolution trend is similar 
to the results presented in Figure 8 for PUC, the time to attain steady state solution 
slightly increases. Evolutions of the selected effective moisture terms is plotted in Fig- 
ures 12(c) (d) confirming again a negligible difference between steady state and transient 
analyses. While the degree of heterogeneity is pronounced only slightly, the influence of 
initial conditions on the prediction of effective properties is significant. 
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Figure 12: Evolution of local moisture (SEPUC): (a) loading conditions I, (b) loading conditions II; 
Evolution of homogenized material parameters as s function of time: (c) ^ ^, (d) C^, 

Multi-scale analysis 

Combining all the previous results suggests that for a reasonably small macroscopic 
time increment (from one to two hours sufficient to reach the steady state conditions on the 
meso-scale) the influence of non-local terms in Eqs. (22) and (23) should be negligible. 
Therefore, the macroscopic response should be invariant with respect to the adopted 
analysis carried out on meso-scale whether the transient or steady state. This is evident 
from the results plotted in Figure 13 showing evolution of macroscopic temperatures and 
relative humidities at selected nodes of mesh, see Figure 5. 




Figure 13: Comparison of different macrostructural computations at selected nodes - (a) temperature 
evolution, (b) moisture evolution 

This result is quite encouraging particularly with reference to the analysis of large 
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historical structures, since avoiding a transient analysis on each mesoscopic unit cell may 
considerably reduce the computational cost. 

4. Conclusions 

The present paper gives a brief overview of two particular aspects of the modeling of 
historical masonry structures, which show a certain degree of irregularity on the meso- 
scale. This observation promotes application of so called Statistically Equivalent Periodic 
Unit Cell being sufficiently simple in comparison to, yet sufficiently representative of, a 
real masonry. This issue was addressed first in Section 2 to see that under steady state and 
linear conditions the degree of heterogeneity, associated with a given meso-structure, may 
not play a significant role in the prediction of effective macroscopic response, recall Table 1. 
On the contrary, the morphological details will become important once considering a 
nonlinear response crucially dependent on the actual distribution of local fields. 

The latter comment was partially examined next in Section 3 devoted to the nonlin- 
ear fully coupled multi-scale analysis of simultaneous heat and moisture transport. The 
principal result of this study is seen in the possibility of deriving the instantaneous, time 
dependent, macroscopic response from a steady-state analysis performed on the lower scale 
reflecting all morphological details. However, keep in mind that this finding is strictly 
valid for the present problem being a collection of the selected climatic conditions, mate- 
rial composition and the assumed constitutive model, and should not be generalized. For 
other cases the theoretically predicted dependence of macroscopic response on the actual 
RVE size may prove non-negligible [8, 23]. 

It is our present interest to exploit the two advancement in computational efficiency 
(SEPUC, steady state meso-scale problem) in the analysis of full scale 3D model of Charles 
Bridge. Special attention will be devoted to the implementation efficiency in the frame- 
work of hybrid parallel computing. 
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